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ABSTRACT 

Aims. In order to maximize the scientific return of the telescopes located at the Paranal Observatory, we analyse the 
properties of the atmosphere above Paranal in the near- infrared (NIR). 

Methods. We estimate atmospheric extinction in the spectral range 1.10-2.30 fim (J, Js, H, and Ks) using standard 
stars that were observed during photometric and clear nights with ISAAC on the Very Large Telescope UTl between 
2000 and 2004. We have built a database consisting of hundreds of observations, which allows us to examine how 
extinction varies with airmass and the column density of water vapour. In addition, we use theoretical models of the 
atmosphere to estimate Rayleigh scattering and molecular absorption in order to quantify their fractional contribution 
to the total extinction in each filter. Finally, we have observed 8 bright red standard stars to evaluate filter color terms. 
Results. We find that extinction coefficients are < 0.1 mag airmass"^ in all the considered bands. The extinction 
coefficient in the J-band strongly depends on the column density of water vapour. Molecular absorption dominates the 
extinction in J, H and Ks , while Rayleigh scattering contributes most to the extinction in Js . We have found negligible 
color terms for J, H and Ks and a non-negligible color term for Js. 

Key words. Atmospheric Effects, Site Testing 



1. Introduction 

The ESO Paranal Observatory is located on the edge of 
the Atacama Coast in Northern Chile, 120 km south of the 
city of Antofagasta, at an altitude of 2635 meters above 
sea leveL The conditions at the observatory are charac- 
terised by a high fraction of clear nights and low levels of 
humidity (Lombardi et al. I2009p . which are important for 
near-infrared (NIR) observations. However, there has been 
little direct measurement of the properties of the atmo- 
sphere above Paranal in the NIR. In order to fill this gap, we 
have calculated the extinction coefficients using photomet- 
ric standard stars observed with the Infrared Spectrometer 
And Array Camera (ISAAC) at the Very Large Telescope 
(VLT) Unit Telescope 1 (UTl) as part of the ISAAC cali- 
bration plan. In addition, we have estimated the fractional 
contribution of Rayleigh scattering, molecular absorption 
and aerosol scattering to the extinction. This estimation 
required the application of a theoretical approach retrieved 
from the literature. Finally, to complete the analysis, we 
also calculated the instrumental color terms using 8 red 
stars observed under photometric conditions. 



2. Database and data reduction 

Our dataset consists of standard star observations covering 
a period of 5 years (from March 2000 through to December 



2004). The observations were made during clear nights with 
the ISAAC short-wavelength (SW) arm in J, Js, H, and 
Ks- All data, including the relevant calibration frames, were 
retrieved from the ESO archive. Clear nights were identified 
by inspecting the nightly observatory weather reports. All 
stars have magnitudes on the LCO system (Persson et al. 
HMD- 
Table [T] lists the characteristics of ISAAC filters. 
(ISAAC User Manual HOMl) , while in Table El at the end 
of this Research Note, we note all the standard stars used 
in the analysis. The data reduction has been performed us- 
ing the ISAAC Pipeline. Each frame has been corrected for 
electronic artifacts, dark subtracted and fiat fielded. Each 
standard star is imaged over a grid of five positions, one just 
above the center of the array and one in each quadrant. The 
pipeline computes a set of instrumental magnitudes which 



Table 1. Properties of ISAAC SW imaging filters used in 
this study. 



Filter 


Central Wavelength [^iin] 


Width [^lm] 


Width [%] 


J 


1.25 


0.29 


23 


Js 


1.24 


0.16 


13 


H 


1.65 


0.30 


18 




2.16 


0.27 


13 



Send offprint requests to: G. Lombardi 
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Fig. 1. Evolution in time of the zeropoints for J, Js, H and Kg. The solid vertical lines indicate Ml recoating events, 
while dotted vertical lines indicate ISAAC interventions. 



are averaged to deliver the zeropoint {ZP) uncorrected for 
extinction. 

ZP uncertainties are computed as cr = azPi/ \/n, where 
n is the number of times, usually five, a single standard 
is observed, and cr^p. is the scatter about the mean ZP 
reported by the pipeline. In our analysis, we have rejected 
points having a > 0.050 mag. Our final sample contains 
575 data points in J, 603 data points in Jg, 604 data points 
in H and 667 data points in Kg. 

3. Data analysis 

3.1. Time evolution of the zeropoint 

In Figure [1] the zeropoints of the four bands are shown as 
four time series. Different colors represent different bands. 
We clearly see that, for each band, the zeropoints are char- 
acterized by trends within well defined time intervals. A 
further inspection demonstrated that the intervals are de- 
limited by technical or maintenance interventions, either on 
the instrument or the telescope (e.g. Ml recoating). This 
means that, in the mentioned intervals (or periods, P), the 
computed zeropoints are affected by deterioration of the 
telescope optics after a recoating of UTl primary mirror 
(Ml) or instruments troubles (ISAAC technical interven- 
tions). A list of the events that occurred between 2000 and 
2005 is reported in Tabled 



The deterioration in the aluminum coating of the tele- 
scope mirrors is mainly due to dust and oxidation. It results 
in a progressive reduction in the mirror reflectivity, an in- 
crease in the thermal background emission (Frogel I1998P 
and a consequent decrease in the zeropoints. A change (in- 
crease or decrease) in the zeropoint can also be due to inter- 
ventions on ISAAC, as they affect the instrument configu- 
ration. Our idea was to consider each period separately and 
to remove, from each one of them, the time dependency by 
subtracting a linear fit of the zeropoint with time. In prac- 
tice we adjust all zeropoints to the value of the intercept 
of the fit at the beginning of the considered period. During 
the analysis we have confirmed that a linear fit was accurate 



Table 2. Technical events occurred between 2000 and 2005. 



Item 


Date 


Event 


PI 


2000 02 10 


Ml recoating 


P2 


2001 03 27 


ISAAC intervention 


P3 


2001 08 21 


ISAAC intervention 


P4 


2001 10 14 


ISAAC intervention 


P5 


2002 03 14 


ISAAC intervention 


P6 


2002 11 25 


Ml reeoating 


P7 


2003 03 28 


ISAAC intervention 


P8 


2004 01 26 


ISAAC intervention 


P9 


2004 04 03 


Ml reeoating 
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Fig. 2. ZPq for each period. 

enough, therefore higher polynomials have not been used. 
The fits for periods P3 and P6 are not well constrained, 
because the small number of points within these two peri- 
ods. Therefore, they are not been considered further in this 
analysis. 

The evolution in the zeropoint within each period is 
fitted with 

ZP{t)^Ct + ZPq (1) 

where t is the time, ZPq corresponds to the zeropoint at 
the beginning {t = 0) of the considered period, and C the 
slope of the fit. The differences ZP{t = 0) — ZPi have been 
computed and added to each observed ZPi. 

Finally, in each period, we subtracted the intercept of 
the fit at i = from the corrected ZPi. In this way we 
have eliminated both the offset between periods and the 
evolution with time within a period. This represents a clear 
advantage. We can now work with a single dataset, rather 
than 9. 

It is interesting to check the variation in ZPq after each 
event. If the event did not affect the instrument perfor- 
mance, ZPq values should be all identical. As shown in 
Figure m we observe significant offsets between events. 

3.2. Effects of the Precipitable Water Vapour 

The transmission of the atmosphere in the NIR is greatly 
affected by the presence of water. Considering the classical 
Johnson filters, an increase in the column density of water 
vapour causes (J — K) and {H — K) to be bluer, because 
the effect of water vapour is stronger in K (AA = 0.41 
/im) than in J and H (Frogel 119981 see also Figure I in 
Manduca&Bell 1979). On the other hand, the ISAAC 
filter has a narrower bandwidth (AA = 0.27 /xm, see Table 
[1]) and avoids the strongest water lines. Therefore we ex- 
pect that, in our case, J and H will suffer more from large 
amounts of water in the atmosphere. 

The amount of water above Paranal has been contin- 
uously monitored since July 2000 using the images from 
the Geostationary Operational Environmental Satellite 



Fig. 3. Monthly trend of the PWV on Paranal in photo- 
metric and clear nights. 



(GOES). The sampling is every 3 hours starting at OOh 
UT. 

To retrieve the Precipitable Water Vapour (PWV) 
from GOES images a model exists developed by A. 
Erasmus under contract with ESO (Erasmus&Pcterson 
1997; Erasmus&SarazinHnnOl Erasmus&Sarazin2002). The 
method is based on combining satellite current image in the 
6.7 /im and 10.7 /zm channel with wind and temperature 
profiles forecasted by a prediction center. 

Typical satellite observations at about 6.5 /im are sen- 
sitive to emissions from water vapour resident in the layer 
between about 600 mbar (~4400 m) and 300 mbar (-9000 
m). For what concerns GOES vertical resolution, the 6.7 
/^m channel is located near the center of a strong water 
vapor absorption band and under clear sky conditions it is 
primarily sensitive to the relative humidity averaged over 
a depth of atmosphere extending from 200 to 500 mbar 
(Soden&Bretherton 1993). The horizontal resolution is suf- 
ficient since the atmosphere surrounding Paranal is not ex- 
pected to vary significantly over an area that extends sev- 
eral lO's of km. 

Figure [3] shows the yearly trend in the amount of PWV 
sampled at Paranal during photometric and clear nights. 
We do see a yearly periodic modulation with very high val- 
ues (up to 13 mm in 2002) in the trimester January-March 
corresponding to the so called Bolivian winter. As shown in 
the figure, the median PWV at Paranal during photometric 
and clear nights is 2.3 mm. 

Figure |4] shows the relationship between zeropoints and 
the PWV in the considered bands. For each band we have 
calculated a weighted linear fit. The fits are repeated twice 
more after rejecting 3-tT outliers. In Table |3] we report the 
slopes and the RMS of the fits. For Js, H , and Kg the trend 
with the column density of PWV is quite slight, while for J 
it is very significant. The scatter about the best fit for J is 
between 2 and 4 times larger that the scattered measured 
for the other filters. 
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Table 4. Calculated k in [mag airmass^^] for J, Jg, H and Kg at Paranal. 



Band 


PWV-corrected ZP 


PWV 


= 0-2 mm 


PWV = 


2—4 mm 


PWV 


= 4—7 mm 




/t 


disp. 


K 


disp. 


K 


disp. 


/t 


disp. 


J 


0.072 


0.040 


0.038 


0.038 


0.060 


0.036 


0.090 


0.033 


Js 


0.048 


0.019 


0.034 


0.019 


0.040 


0.020 


0.058 


0.018 


H 


0.034 


0.015 


0.035 


0.015 


0.030 


0.015 


0.053 


0.016 


Ks 


0.043 


0.013 


0.040 


0.015 


0.046 


0.014 


0.042 


0.012 




PWV [mm] 

Fig. 4. Weighted linear fit of the zeropoints as a function 
of the PWV in the four bands. 

Table 3. Zeropoints- versus-PWV slope and dispersion 
(RMS). 



Band 


Slope 


Dispersion 




[mag per unit of PWV] 


[mag] 


,7 


-0.007 ± 0.001 


0.048 


Js 


-0.003 ± 0.001 


0.024 


H 


-0.003 ± 0.001 


0.019 


Ks 


-0.002 ± 0.001 


0.011 



4. Determination of the atmospheric extinction 
coefficients 

The relationship between ZP and airmass {X) defines the 
extinction curve (Bouguer curve). The extinction coefficient 
(k) is the slope of a linear fit to this curve. The fit is done 
three times. Between each fit, 3-cr outliers are rejected. 

The extinction coefficients have been calculated in two 
different ways: 

a) correcting the zeropoints for the amount of PWV before 
computing k; 

b) considering separately the zeropoints in different PWV 
ranges in order to evaluate function of the PWV. 

In case a we offset the zeropoints using the fits calcu- 
lated in Sj321 and reported in Figure |4] and Table |3l The 
offset to each zeropoint is 

/^ZP,^ ZP,{PWV ^2.Z)~ ZP,{PWV) (2) 




X 



Fig. 5. Bouguer curves computed for the PWV rescaled 
zeropoints in Js, H and Ks- 

where the index i refers to one of the four filters. The ze- 
ropoint is corrected to the median amount of PWV. The 
PWV-corrected zeropoint, rZPi, is then 

rZP, = ZP, + AZP, (3) 

In Figure [5] we have plotted the weighted linear fits to the 
curves for the rZP in the four bands. The extinction coef- 
ficients are reported in column 2 of Tabled The dispersion 
(RIVIS) of the points around the best fits are also reported. 

In case b we have calculated the extinction coefficients 
for three different PWV ranges: PWV = 0-2 mm, PWV = 
2-4 mm, and PWV = 4-7 mm. The results are reported 
in columns 3, 4, and 5 of Table S] and plotted in Figure 
El As expected, the J-band extinction coefficient has the 
greatest sensitivity to the amount of PWV. It increases by 
~0.05 mag airmass"^ over the range of PWVs considered 
here. The increase in the extinction coefficients for Js and 
H are less (~0.03 and ^0.02 mag airmass"^, respectively), 
while there is a negligible effect in Kg- 

The dependence of the extinction coefficients on the 
amount of PWV underlines the importance of observing 
IR standards (more important for J and less important for 
Ks) at roughly the same airmass and at roughly the same 
time as the science target. 

5. The fractional contribution to the extinction 

According to Hayes&Latham (|1975p the largest contrib- 
utors to atmospheric extinction are Rayleigh scattering 
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Table 5. Raylcigh scattering coefRcients and vertical 
Rayleigh extinction in [mag airmass""'^] calculated for 
Par anal. 





B(T, P, A) 




J 


0.106 ± 0.011 


0.043 ± 0.005 


Js 


0.109 ± 0.011 


0.046 ± 0.005 


H 


0.035 ± 0.004 


0.005 ± 0.001 


Ks 


0.012 ± 0.001 


0.001 ± 0.001 



Table 6. Calculated vertical molecular extinction in [mag 
airmass^^l for Paranal. 



Molecule 


.7 


Js 


H 




CH4 


0.001 


0.000 


0.004 


0.015 


CO2 


0.000 


0.001 


0.005 


0.040 


H2O 


0.189 


0.018 


0.027 


0.004 


N2O 


0.000 


0.000 


0.000 


0.001 


O2 


0.004 


0.007 


0.000 


0.000 




0.194 


0.026 


0.036 


0.060 



Fig. 6. Extinction coefficients in J, Js, H and Kg calcu- 
lated for PWV = 0-2 mm, PWV = 2-4 mm, PWV = 4-7 
mm. The dotted line represents the k calculated with the 
PWV rescaled ZP, while the gray zone is the dispersion 
(see column 2 of Table S]). 

(K-nay), molecular absorption {k„ioi) and aerosol scatter- 
ing (Kaer)- The efhcicncy of aerosol scattering from parti- 
cles that are smaller than a few microns at A > 1.0 /im is 
negligible (Lombardi et al. I2008p . therefore for the ISAAC 
filters (Tabled)) we have computed the vertical atmospheric 
extinction coefficients as the sum of KRay and Kmoi 

K\ = KRay.X + f^mol.X (4) 

Variations in KRay and Kmo/ with time are linked to oc- 
casional and periodic climatic changes in the atmosphere 
above the site (Lombardi et al. I2009P . For this reason, the- 
oretical models are valid only to calculate average contri- 
butions to the extinction coefficients. 

5.1. Rayleigh scattering 

Rayleigh scattering of unpolarised light is due to particles 
with dimensions that are ^ A. It has been extensively dis- 
cussed by Penndorf p957p . According to Hayes&Latham 
(fW5)) . Rufener ((TO5S|) and Burki et al. (ITM5|) the vertical 
Rayleigh extinction can be expressed as 

KRay,\ = BX^^ (5) 

where the Rayleigh scattering coefficient is a complex 
function of the refractive index, n{T, P, A), and the number 
of particles, N{T,P). T and P are the site's mean tem- 
perature and pressure. B can be written as B{T,P,X) and 
expressions for n, N and B are given by Penndorf (1957). 

In the years between 2000 and 2004 a mean tempera- 
ture of (12.8± 0.5)°C and a mean pressure of (743.5 ±02) 
hPa have been measured at Paranal (Lombardi et al. I2009P . 
We have used the expressions in Penndorf p957p and 
Hayes&Latham (|1975p to compute B{T,P,X) and KRay,\ 
for the four ISAAC filters. The values are reported in Table 

El 



Table 7. Vertical atmospheric extinction coefficients [in 
mag airmass"^] and fractional contributions of Rayleigh 
scattering and molecular absorption. 





ka 


fraction of Kjiay,x 


fraction of Kjy^ol,x 


J 


0.237 ± 0.010 


18% 


82% 


Js 


0.072 ± 0.010 


64% 


36% 


H 


0.041 ± 0.010 


12% 


88% 


Ks 


0.061 ± 0.010 


2% 


98% 



As expected, Rayleigh scattering is more efficient in J 
and Js. It is an order of magnitude less efficient in H and 

5.2. Molecular absorption 

The contribution to the extinction originating from molec- 
ular absorption was modeled using an IDL driver to the 
Reference Forward Model (RFM). RFM is a GENLN2- 
based line-by-line radiative transfer code developed by Ann 
Dudhia at the Atmospheric, Oceanic and Planetary Physics 
Institute at Oxford University (UK) to analyse data from 
MIPAS on-board ENVISAtQ. The code was run using the 
2008 version of the HITRAN database (Rothman et al. 

mm- 

For the current work, only lines caused by H2O, CO2, 
N2O, CH4 and O2 were considered. We slightly modified 
the tropical atmospheric profile available with the RFM 
code to carry-out simulations corresponding to a temper- 
ature of 12.8°C and an atmospheric pressure of 743.5 hPa 
(Lombardi et al. I2009p . The amount of water vapour was 
scaled to match the median PWV for Paranal (2.3mm). 
The code also allows one to vary the airmass of the line-of- 
sight or to select a specific molecule. The extinction in each 
filter was estimated by multiplying the atmospheric trans- 
mission produced by the model with the filter transmission 
curve (see ISAAC User Manual r2009p . Results are reported 
in Table [61 



http:/ /www.atm.ox.ax.uk/RFM 
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Table 8. Observed Persson's red standard stars for the determination of the color terms (the computed Cm is also 
reported). 





R.A. 


Docl. 












Star 


(J2000) 


(J2000) 


,7 - 


J - Jis 


./ - Js,is 


H - His 


Ks - Ks.is 


cskd-21 


12 32 10.9 


-63 43 16 


0.724 


-0.262 ± 0.022 


0.018 ± 0.022 


0.192 ± 0.024 


0.722 ± 0.032 


cskf-12 


12 31 30.1 


-63 51 03 


1.110 


-0.281 ± 0.017 


-0.011 ± 0.015 


0.204 ± 0.011 


0.731 ± 0.012 


cskf-14a 


12 31 45.9 


-63 49 36 


1.734 


-0.281 ± 0.023 


-0.005 ± 0.019 


0.207 ± 0.018 


0.708 ± 0.022 


cskd-20 


12 32 04.0 


-63 43 46 


1.826 


-0.284 ± 0.026 


0.025 ± 0.026 


0.219 ± 0.017 


0.725 ± 0.020 


cskd-9 


12 31 16.7 


-63 40 11 


2.219 


-0.243 ± 0.039 


0.047 ± 0.032 


0.187 ± 0.030 


0.690 ± 0.033 


cske-23 


12 31 56.0 


-63 37 43 


2.735 


-0.244 ± 0.027 


0.026 ± 0.025 


0.221 ± 0.017 


0.736 ± 0.020 


cskd-34 


12 31 23.6 


-63 46 45 


2.826 


-0.286 ± 0.021 


0.034 ± 0.022 


0.198 ± 0.019 


0.725 ± 0.023 


cskd-16 


12 31 57.8 


-63 42 21 


2.862 


-0.254 ± 0.021 


0.056 ± 0.020 


0.203 ± 0.019 


0.732 ± 0.021 








C„ 


0.008 ± 0.010 


0.020 ± 0.009 


0.003 ± 0.006 


0.002 ± 0.008 



5.3. Fractional contributions 

In column 1 of Table [71 we list kx for J, J, , H and Kg . 
The values have been derived using equation @ and the 
computed values for Kiiay,\ and Kmoi.x- The uncertainties 
derive from the propagation of the uncertainties of T and 
P in Lombardi et al. ( 2009[ ) and as to be assumed has an 
upper limit to the variation of the theoretical k\. 

For Js, H and Kg, there is good agreement between 
the theoretical k\ and the empirical k calculated from the 
PWV rescaled zeropoints (column 2 of Table H]). For the 
J-band, there is a significant discrepancy, which is due to 
the high opacity of the atmosphere at the red end of the J 
filter. The red end of the J filter is effectively defined by the 
atmosphere. Changes in the amount of PWV in the atmo- 
sphere shifts this edge, which leads to larger photometric 
uncertainties, as demonstrated by the large scatter about 
the best fits in Figures S] and [S] As noted in 21 this result 
demonstrates the importance of observing IR standards at 
roughly the same airmass and at roughly the same time as 
the science target when aiming for precise photometry in 
the J-band. 

Columns 2 and 3 of Table [7] show the fractional contribu- 
tions of KRay,x and Kmoi.x to the total vertical atmospheric 
extinction. Molecular absorption dominates the absorption 
in J (82%), H (88%) and (98%). Molecular absorption 
and Rayleigh scattering contribute almost equally to the 
absorption observed in Jg. This is due to the design of the 
ISAAC Js filter which is narrower than J and avoids the 
strong atmospheric absorption lines at the end of J-band. 

6. Determination of the color terms 



According to Amico et al. 



, the J, H and Kg filters 



of ISAAC closely match the filters tabulated in Persson et 
al. (jl998p . We therefore expect negligible color terms for 
these filters. The Js filter of ISAAC, on the other hand, is 
significantly different to J, so we expect a significant color 
term for Jg. 

We observed 8 bright red stars from Persson et al. (|1998p 
with ISAAC in J, Js, H and Ks under photometric condi- 
tions (see Table HI). Each star has been observed over a grid 
of four positions (one for each quadrant of the detector) 
with a windowed detector and short integrations (< Is) 
to avoid saturation. The integrations were shorter than the 
minimum integration time required for full readout of the 
array, so only a subarray of the detector was read out. 
We integrated the flux in apertures of fixed size to derive 
instrumental magnitudes for each star. The final instrumen- 
tal magnitude associated to the star (mis) is the weighted 



-0.1 
■ -0.2 
-0.3 



0.1 



-0.1 - 



0.3 



0,1 
I 0.8 



0.7 - 
0.6 - 







(J - Ks) 



Fig. 7. Plots of {rricat — 'm'is)i versus (J — Ks)i. The slope 
of the linear fit is Cm and is tabulated in Table HI 



average of the four instrumental magnitudes corrected us- 
ing the atmospheric extinction coefficients determined in Sj4] 
(see column 2 in Table H]). 

In Figure[7]we plot {Mcat — 'niis)i against {J—Ks)i and 
fit for the slope Cm- An arbitrary constant vertical offset 
has been applied to the points in each of the graphs. 

Table[5]lists {meat — mis)i for J, Js, H and Ks and the 
color (J — Ks)i. The slopes. Cm, are also reported. 

As expected, Cm does not differ from for J, H oi Ks, 
while Cm = 0.020 ± 0.009 for Js. 



7. Conclusions 

In this Research Note, we have characterized the extinction 
properties in the NIR at the Paranal Observatory (Chile). 
We have used standard stars observed with ISAAC in J, 
Js, H, and Ks during photometric and clear nights at the 
Very Large Telescope UTl between 2000 and 2004. 

For each star we calculated the zeropoint uncorrected 
for extinction using the ISAAC Pipeline. A correction is 
then performed on the whole dataset in order to eliminate 
the affects on the zeropoints due to technical events (such as 
recoatings of the primary mirror and ISAAC interventions). 
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We used two different methods to derive extinction co- 
efficients. In the first method, we account for the varition 
in the ZPs with the PWV by rescaling the zeropoints to the 
median PWV (2.3 mm) measured at Paranal between 2001 
and 2005. We then fit a linear relation to Bouguer curves 
to determine the extiction co-efficients. We have obtained 
Kj = 0.072±0.040, Kj^ = 0.048±0.019, kh = 0.034±0.015 
and kk^ = 0.043 ± 0.013 mag airmass"^. In the second 
method, we have calculated k considering the zeropoints 
in three different PWV ranges: PWV = 0-2 mm, PWV = 
2-4 mm, and PWV = 4-7 mm. As expected, the extinc- 
tion coefficient in J is more sensitive to the amount PWV 
than the extinction coefficients of other filters, increasing by 
~0.05 mag airmass"^ between lower and upper ranges for 
the PWV. For comparison the coefficients for Jg and H in- 
crease by ^-^0.03 mag airmass"^ and ~0.02 mag airmass"^, 
respectively, while there is negligible change for Kg- 

Using a theoretical approach, we have found that molec- 
ular absorption contributes most to the total absorption in 
J (82%), H (88%) and (98%), while Rayleigh scattering 
contributes most to the total absorption in the J^-band. 

We have calculated the color terms using 8 bright stars 
observed in photometric conditions. We have found negli- 
gible color terms in J, H and Kg and a non-negligible color 
term in J, (0.020 ± 0.009). 
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Table 9. Observed Persson's standard stars from the ISAAC Calibration Plan. 
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